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Abstract 

We investigate the spin structure function g2{x, Q^) in the framework of 
the operator product expansion and the renormahzation group. The twist-3 
operators appearing in QCD are examined and their relations are studied. 
It is noted that operators proportional to equation of motion appear in the 
operator mixing through renormahzation, which can be studied from the 
relevant Green's functions. We also note that the coefficient functions can 
be properly fixed after the choice of independent operators. 
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Recently there has been much interest in the nucleon spin structure functions which 
can be measured by the polarized deep inelastic leptoproduction. The nucleon spin 
structure is described by the two spin structure functions gi{x, Q^) and g2{x, Q"^). The 
data from the EMC Collaboration on gi have prompted many authors to reanalyze 
the QCD effects on this process and recent data from SMC at CERN p| and E142 
at SLAC have attracted much more attention in connection with the Bjorken sum 
rule. On the other hand, the experiment on g2 is expected to be performed at CERN, 
SLAC and DESY in the near future. 

For gi, the leading contribution comes only from twist-2 operators in the operator 
product expansion (OPE), while for g2 the twist-3 operators also contribute in the 
leading order of l/Q"^ in addition to the twist-2 operators 0]. This leads to new 
features which do not appear in the analyses of other structure functions. 

After the old papers ^, |^ which discussed the QCD effects on the polarized 
process, the above problems have been addressed by many authors It turns out 

that the procedure to obtain the QCD corrections to the Wilson's coefficient functions 
of each operator is not so straightforward for g2. The presence of the operators pro- 
portional to the equation of motion brings about some complication for the twist-3 
operators, through the course of renormalization, which is a general feature character- 



istic to the higher-twist operators ||10|. 

In this paper we focus our attention on the structure function g2{x, Q"^), and study 
twist-3 operators in QCD in the framework of OPE and renormalization group. We ex- 
amine the operator mixing under the renormalization and point out that the coefficient 
functions depend upon the choice of the basis of the independent operators. 

Let us start with the basic formalism based on the operator product expansion for 
analyzing the polarized leptoproduction The hadronic tensor Wf^i, is given by the 
absorptive part of the forward virtual Compton amplitude. It's anti-symmetric part 
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can be expressed by the two spin-dependent structure functions gi and g2 as : 



where q is the virtual photon momentum and p is the nucleon momentum, x is the 
Bjorken variable x = Q'^/2p ■ q = Q"^ /2Mv , p ■ q = Mv and q^ = —Q^. M is the 
nucleon mass and = u{p, s)'y^'y5u{p, s) is the covariant spin vector. 

Applying OPE to the product of two electromagnetic currents, we get for the anti- 
symmetric part, 

j 



where i?"'s are the composite operators and i?"'s the corresponding coefficient func- 
tions. In (|l|), Rq represents the twist-2 operators and the other operators inside the 
summation over j are the twist-3 operators. For simplicity, let us consider the flavor 
non-singlet case. (In the following expressions, we suppress the flavor matrices Aj for 
the quark field ) The explicit forms of the twist-2 operators are given by 

^a/xi-Mn-i ^ ^'^-i^^g^^/^Mi . . . Dt^n-i}^ _ (traces) , 

where { } means symmetrization over the Lorentz indices and — (traces) stands 
for the subtraction of the trace terms to make the operators traceless, which will be 
suppressed in the following. For the twist-3 operators, we have 



Zip — 



n 



(n - l)ip-f5Y D^'"' ■ ■ ■ Di'^-^^ip 

n-l_ 

- E i^lbl^'D^^D^' ■ ■ ■ D'"-W'+^ ■ ■ ■ L)^"-i^V^] , (2) 
1=1 

R^l■■■^^n-l = z'^-2rri^757'^DiA*i . . . L)''"-27'^"-i}7/>, (3) 

Rr"^"~' = ^{Vk-Vn-l-k + Uk + Un-l-k), (4) 

2n 



where m in represents the quark mass (matrix). The operators in contain the 
gluon field strength G^j, and the dual tensor G^^ = ^e^^afsG"^ exphcitly and given by 

Uk = e-^gSi^D'"' ■ ■ ■ G"''^ ■ ■ ■ D'^'^-^j''"-'^, 

where S means the symmetrization over fii and g is the QCD coupling constant. 

It is by now well-known p, Q that the operators (@ - are related through the 
operators which are proportional to the equation of motion ( EOM operators ), 

r} — 1 

It is not difficult to obtain the following relation : 

r? - 1 

^aMi -/^n-i ^ >^ + y (n - 1 - k)Rl^'-^"-' + ••'^"-^ (5) 

This equation means that these operators are not all independent but constrained by 
this relation, and therefore the renormalization procedure for the twist-3 operators 
becomes rather complicated. Although we realize that physical matrix elements of 
the EOM operators vanish |jTO|, we keep them to study the operator mixing through 



renormalization because their off-shell Green's function do not vanish and are relevant 
for our analyses shown below. 

The anomalous dimensions which enter into the renormalization group equation 
for the coefficient functions, are obtained from the renormalization constant for the 
composite operators. The anomalous dimensions of twist-3 operators in QCD at the 
one-loop level have been studied in ref. Here we reanalyze the operator mixing 
problems by keeping the EOM operators. We examine what happens to the renor- 
malization if there are several operators which are related by constraints. Let us first 
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illustrate the points with the scalar field theory for simplicity. Assume that we have 
two composite operators Ri and R2 which are related through the equation of motion: 

Ri = R2 + E, (6) 

where £" is a EOM operator. The general form of E will he E = v4(</)) 5S/5(j) where S is 
the action and A{(f)) is a some function of 0. Since we have three composite operators 
with a constraint @, we can choose any two operators among three as an independent 
basis. 

First we choose Ri and R2- Then renormalization takes the following form: 



R2 1^ ^21 ^22 J \R2 



(7) 



Here the suffix R (B) denotes the renormalized (bare) operator. Next, we change the 
operator basis to Ri {i = 1 or 2) and E. Then the renormalization matrix becomes 
triangular [jlT], 

(8) 



\ = ( \ ( Ri \ 

E In \ ^ Zee)[Re), 



since the counterterm to the operator E should vanish by the equation of motion. 
From (|^ and (^), we can derive the following relations between the renormalization 
constants Z's: 

= Zii + Z12 = Z21 + Z22 = Z22, (9) 
Zee = Z22 — Z\2 = ZxY — Z21 , Zie = —Zyi , Z2E = Z2i- 

Eq.(^ is a consistency condition for the renormalization of the operators satisfying (|]). 

Here we note that the arbitrariness in the choice of the operator basis does not 
enter into the physical quantities. What we really need is the physical matrix element 
of operators. This fact can be explicitly confirmed by noting that the physical matrix 
element of the operator E vanishes ||10| , |ll[: {p\E\p') = 0. Taking the physical matrix 



element of we conclude that the renormalization constant is only relevant. We 
reach the same conclusion also from (|^ by taking account of the fact: {p\Ri\p') = 
{p\R2\p') and @. The explicit one loop calculation shows that for the simplest twist-3 
operator ipD^Dj^tp, as discussed by Jaffe P], the above argument actually holds [|T^ . 

Now we go back to the case of polarized deep inelastic scattering and give explicit 
results at the one-loop level for n = 3 as an example. In this case we have four operators 
with the constraint (15|). 



Rf 



\Rm + Ri + R, 



eqi 



(10) 



where Lorentz indices of operators are suppressed. Now we choose Rp, Rm, and-Ri as 
independent operators and eliminate EOM operator R^q. In order to renormalize above 
operators, it turns out that a gauge non-invariant EOM operator R^qi comes into play, 



Its presence is allowed because it vanishes by the equation of motion [rT|. We found 
the following renormalization constant for the composite operators. 



(11) 
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where Zij are given in the dimensional regularization D = A — 2e: 



Zij — Sij 
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e IGtt 



with 



iC2{R) + fC2{G), 



\C2{R) + fC2iG), 



(12) 



— |C2(-R) + |C2(G), Z12 '- 

z,3 = 3C2{R) - lC2{G), Zu ^ 
Z21 = qC2{R) — gC2(G'), ^;22 - 

Z23 = —^C2{R) + Y2^2(G), ^24 = -gC2{G)^ 

Z33 = 6C2{R), Zu = 0. 

The quadratic Casimir operators are G2{R) = 4/3 and C2{G) = 3 for the case of 

QCD. In the above calculations, it is necessary to compute the off-shell Green's func- 



tions of the composite operators. Otherwise some informations on the renormahzation 
constants associated with the EOM operator will be lost. 

As explained before, one can choose other operators as an independent basis using 
([T0|) . The consistency conditions corresponding to @ read in this case, 

Zu + Zi2 = Z21 + Z22, ^Zii + Zi3 = 1^21 + ^23 + |^33; 

^13 ~ 3^12 = Z23 — -^Z22 + 3^33- 

These equalities are indeed satisfied with ([I^). 

The authors in ref.p] have discarded the fermion bilinear operator j^^/^i ' ^*"-! 
their analyses. For the present case of n = 3, the renormahzation constant matrix 
becomes 
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where Zij are defined in (0). In this basis, our results agree with those in ref. P]. 
Although this choice of basis might be economical especially for general n, it is never 
compulsory. One can eliminate any operator as well. Any choice of basis leads to a 
unique prediction for the moment with the properly determined coefficient functions 
discussed below. 

To make a prediction for the moments, we must determine the coefficient functions 
corresponding to the "hard part" of the process. The coefficient functions at the tree 
level can be obtained by considering the short distance expansion of the current product 
in the presence of external gauge field A^. Up to twist-3 operators, the OPE reads p|, ^ 

I J d^xe*'?"T(J^(a;)J,(0))^ 

n=l,3,-- \V / 

Here if we use (^, i?^^^ can be eliminated in terms of other operators. From 
this observation, it is inferred that the (tree-level) coefficient functions depend upon 

6 



the choice of the independent operators. In the basis of independent operators which 



includes Rp, we conclude that at the tree- level, 



K(tree) = E^(tree) = 1 , Ell^{tTee) = ^"(tree) = . 



(15) 



On the other hand, if we ehminate Rp we have, 

n — 1 

i?,"(tree) = l, i?;;(tree) = 



n 



E"(tree) =n-l-j. 



(16) 



Now we shall write down the moment sum rule for g2. Define the matrix elements 
of composite operators between nucleon states with momentum p and spin s by 



{p,s\RT-''"-'\p,s 
{p,s\RT-^--'\p,s 



{p,s\Rr-^-^\p,s 



{p,s\r:i;^-'^-'\p,s 



n — 1 



n 



0. 



(17) 
(18) 
(19) 
(20) 
(21) 



Our normalization in the above definition is such that for free quark target a„ = (i„ = 
Cn = 1. On the other hand, = 0{g^). Using (|l^ - ^T]), we can write down the 
moment sum rule for g2. 



Jo 



with the following constraint from 



n — 1 r 



2n 



+ - 



n-2 

enEl{Q') + Y.fiE^{Q' 



(22) 



n — 1 



n 



n — 1 



n 



n-2 



It should be stressed that the explicit form of the evolution of each term in (p2|) 
depends on the choice of operator basis. The anomalous dimensions as well as the 



coefficient functions take different forms depending on the basis. However the 
moment itself remain the same as one expects. For the case of n = 3, we have exphcitly 
checked that any choice of the basis leads to the same expression for the moment. 

Here we note that the coefficient functions for the EOM operators do not mix with 
those for the other operators. This is due to the triangular structure of the anomalous 
dimension matrix in their renormalization-group equation. This fact is essential to the 
equivalence of the moments for any basis. 

Let us now consider the moments of g2 in the basis of Rj and Rm- In this case, 

T7 — 1 1 r "^"^ 1 

Mn = -^a^E-{Q') + 2 hnEHQ') + ^ fiE^iQ^)] . (23) 

with the coefficient functions given in (]T6|). We shall show below that ( p3D indeed holds 
in the leading order of InQ^. 

The operator mixing of gluon-field dependent operators Rj and the mass dependent 
operator Rm through the renormalization reads neglecting the EOM operators, 

y Rm Jj^ y Zmm J y Rm j 

By evaluating the Green's functions of these operators with incoming and outgoing off- 
shell quark states, we obtain the off-diagonal element of the renormalization constant 
matrix Zj^ from which the anomalous dimension reads. 



^ \ = ^ . (24) 



This result is in disagreement with the one given in the fifth reference in p|. 

Now note that the Compton scattering amplitude off the on-shell massive quark 
target has been calculated 0, in the leading order of InQ^, to be 
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Whereas, using the perturbative solution for the renormahzation-group equation, the 
right-hand side of ( p3D becomes for the quark matrix elements, 

RHS of (23) = (^(-^)7;£^."(tree) + ^7^,^.(tree)) InQ^ + ■ ■ ■ . (26) 

Putting 7^ of ( p4|) and the following expression for the anomalous dimensions. 



r\ ™ 1 1 1 

7? = 2C,(i?) [l - ^^^^ + -] , llm = 8C,{R) E -. 



into the above equation with the tree level coefficient functions (|T^), we find (^ 



coincides with (|25|) . Thus we confirm that our results (|24| , p!6D are consistent with the 
moment sum rule. 

Finally let us comment on the Burkhardt-Cottingham sum rule which corresponds 
to the first moment of g2. We see from (0-^ that twist-3 operators can not be defined 
for n = 1. Therefore the OPE analysis suggests that the Burkhardt-Cottingham sum 
rule: 

/ dx g2(x, Q"^) = 
Jo 

does not receive any corrections in QCD |T3|, [14 . 



To summarize, we have examined the twist-3 operators in QCD which contribute to 
g2{x, Q^). In the renormalization, there appear operators proportional to the equation 
of motion, which can be studied by computing off-shell Green's functions. Because of 
the relationship among the twist-3 operators, we have to choose a basis of independent 
operators, by which the coefficient functions are properly fixed. We have also noted 
that the Burkhardt-Cottingham sum rule holds in the higher-order of QCD [p^] . 

We expect that future measurements on g2 will clarify the effect of twist-3 operators 
which share some common feature with more general higher-twist effects in QCD. 
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